The spike-specific TCRβ repertoire shows distinct features in unvaccinated or vaccinated patients with SARS-CoV-2 infection

Background The evolving variants of SARS-CoV-2 may escape immunity from prior infections or vaccinations. It’s vital to understand how immunity adapts to these changes. Both infection and mRNA vaccination induce T cells that target the Spike protein. These T cells can recognize multiple variants, such as Delta and Omicron, even if neutralizing antibodies are weakened. However, the degree of recognition can vary among people, affecting vaccine efficacy. Previous studies demonstrated the capability of T-cell receptor (TCR) repertoire analysis to identify conserved and immunodominant peptides with cross-reactive potential among variant of concerns. However, there is a need to extend the analysis of the TCR repertoire to different clinical scenarios. The aim of this study was to examine the Spike-specific TCR repertoire profiles in natural infections and those with combined natural and vaccine immunity. Methods A T-cell enrichment approach and bioinformatic tools were used to investigate the Spike-specific TCRβ repertoire in peripheral blood mononuclear cells of previously vaccinated (n = 8) or unvaccinated (n = 6) COVID-19 patients. Results Diversity and clonality of the TCRβ repertoire showed no significant differences between vaccinated and unvaccinated groups. When comparing the TCRβ data to public databases, 692 unique TCRβ sequences linked to S epitopes were found in the vaccinated group and 670 in the unvaccinated group. TCRβ clonotypes related to spike regions S135-177, S264-276, S319-350, and S448-472 appear notably more prevalent in the vaccinated group. In contrast, the S673-699 epitope, believed to have super antigenic properties, is observed more frequently in the unvaccinated group. In-silico analyses suggest that mutations in epitopes, relative to the main SARS-CoV-2 variants of concern, don’t hinder their cross-reactive recognition by associated TCRβ clonotypes. Conclusions Our findings reveal distinct TCRβ signatures in vaccinated and unvaccinated individuals with COVID-19. These differences might be associated with disease severity and could influence clinical outcomes. Trial registration: FESR/FSE 2014–2020 DDRC n. 585, Action 10.5.12, noCOVID19@UMG. Supplementary Information The online version contains supplementary material available at 10.1186/s12967-024-04852-1.


Background
The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic requires continuous attention focused on the epidemiological, virological, and immunological characteristics of the evolving variants of concern (VoCs).Distinct antigenic characteristics of VoCs could allow the virus to escape from immunity generated through previous infection and/or vaccination, eluding the protective immunity against re-infections and severe disease courses [1].For this reason, it is important to evaluate as many variables as possible that can influence the adaptive immunity against the infection, as well it is relevant to understand how these variables change in relation to the evolution of the epidemiological, virological, and immunological landscape.
SARS-CoV-2 infection and mRNA vaccination were shown to induce spike (S) specific T cells that can recognize and eliminate infected cells [2][3][4][5].These S specific T cells can largely, albeit not completely, tolerate the amino acid mutations that characterize the different VoCs, including Delta and Omicron [6,7], and respond to the viral variants.Such T cell cross-reactivity against VoCs mutant peptides contributes to protection against severe disease, hospitalization, and death, even if the neutralizing antibody response is partially compromised [8][9][10].However, it is important to note that the degree of T cell cross-reactivity may vary among individuals, and the overall effectiveness of current vaccines against new emerging VoCs may be reduced.
The specificity of T cell response is determined by the T cell receptor (TCR), which is produced through a stochastic process of somatic recombination that combines the unique V, D and J gene segments of the TCRα and TCRβ genes (D segments only for TCRβ), resulting in huge range of TCRs with an incredibly diverse antigenic specificity.The third complementarity-determining regions (CDR3s) localized within the TCRα and TCRβ chains are the most hypervariable regions and takes part to direct peptide recognition.Recent advances in TCR sequencing technologies and bioinformatic analysis allow the characterization of the TCR repertoire, that is the collection of diverse and unique TCRs within an individual's immune system, with great throughput and efficiency [11].Analysis of the TCR repertoire provides summary indices of the diversity and clonality of T-cell responses that may be associated with the clinical evolution of a disease, as well as it allows extensive profiling of T cell specificities, despite the complexity of these responses across individuals and groups [11][12][13].
Recent studies found that the TCR repertoire of SARS-CoV-2 specific T cells is highly diverse, which is important for recognizing and clearing the virus [2,14].Moreover, it was observed that COVID-19 patients affected by severe disease have a restricted TCR repertoire and an increased frequency of public TCRs, suggesting that TCR diversity may play a role in determining disease outcomes.The mRNA vaccines elicited a diverse TCR repertoire, indicating a robust T cell response and supporting the efficacy of SARS-CoV-2 vaccines in generating an adaptive immune response [10].These analyses also identified immunodominant TCRs associated with S-specific CD8 + T cell responses [10,15].While these studies demonstrated the utility of T-cell repertoire analysis in identifying conserved and immunodominant peptides with cross-reactive potential among VoCs, predicting disease severity, and informing treatment strategies, they also emphasized the need to extend TCR repertoire analysis to different clinical scenarios.
In this study, we analyzed the TCR repertoires generated during natural SARS-CoV-2 infection in unvaccinated patients or hybrid immunity (infection-induced and vaccine-induced immunity), highlighting distinctive S-specific TCR profiles between the groups.

Study design
An observational longitudinal study was performed on 14 consecutive patients tested positive for SARS-CoV-2 with mild to severe COVID-1912 from January 24th to July 7th, 2022.The characteristics of the participants are summarized in Table 1.The study was conducted according to the standards of the Declaration of Helsinki revised in 2008 (World Medical, 2013), and was approved by the ethical committee of the Calabria Region (Protocol Reference: FESR/FSE 2014-2020 DDRC n. 585, Action 10.5.12, noCOVID19@UMG).Written informed consent was obtained from all the participants before moAbs administration and blood samples collection for the purpose of this study.

PBMCs purification and in-vitro T-cell expansion
Peripheral venous blood was collected in EDTA vacutainer tubes, and peripheral blood mononuclear cells (PBMCs) were isolated by density gradient isolation using Ficoll-Paque (Merck, KGaA, Darmstadt, Germany), according to the manufacturer's instructions.The isolated PBMCs were immediately divided into two aliquots, one subjected to RNA extraction for subsequent TCR sequencing analysis, the other used for the T cell expansion procedure.For in-vitro T-cell expansion, 10 6 PBMCs were seeded cultured with Advanced RPMI Medium 1640 supplemented with 2% human serum, 2 mM L-Glutamine, and 100 U/penicillin/streptomycin in 24-well.Then cells were stimulated for 12 days with 2 μg/mL of SARS-CoV-2 (S-pool) and 50 U/mL rIL-2 and incubated in a humidified CO 2 incubator at 37 ℃ changing medium with fresh S-pool and rIL-2 every 2 days [16,17].The S-pool consisted in 15-mer peptides that overlapped by 10 amino acids and spanned the entire protein sequence of the S protein of SARS-CoV-2 (Uniprot_ID = P0DTC2) [17].

TCR sequencing
TCRβ libraries for NGS sequencing were prepared using the Oncomine ™ TCR Beta-LR Assay (ThermoFisher), according to the manufacturer's protocol.In detail, RNA from PBMCs samples was isolated using the Purelink RNA Mini kit (Thermo Fisher Scientific, Milan, Italy), reverse transcripted through the InvitrogenTM Super-Script IV VILO Master Mix (ThermoFisher Scientific).Libraries were prepared using the TCR beta-LR Assay Kit (ThermoFisher Scientific), which consists of Multiplex AmpliSeq primers target the framework region 1 (FR1) and costant (C) regions of the TCRβ producing a 330 bp 2) with the built-in preset pipeline "Oncomine ™ TCR Beta-LR Assay" [18].

TCR repertoire analysis
TCR repertoire analysis was mostly performed using the Immunarch R package [19].The diversity of TCR repertoires was evaluated by the Gini, Gini-Simpson and d50 diversity coefficients.The Gini coefficient measures the inequality in the frequency distribution of clonotypes, with values close to zero expressing full equality of clonotype frequencies, while a Gini coefficient of 1 reflects maximum inequality between clonotype frequencies, such as the co-presence of hyper-expanded clonotypes and rare clonotypes.The d50 coefficient calculates the minimum number of distinct clonotypes amounting to greater than or equal to 50 percent of a total of sequencing reads obtained following amplification and sequencing.The Gini-Simpson index is the probability of interspecific encounter, i.e., probability that two entities represent different types.Top10, rare and hyperexpanded clonotypes abundance were calculated through the repClonality function of Immunarch.The distribution of Vβ gene segments between COVID-19 groups was performed using the Gene Usage Analysis tool of Immunarch.

In-silico analysis
The GLIPH2 algorithm of turboGLIPH R package was used for clustering of TCRβ sequences [12].The minimum cluster size parameter was set to 8. The mapping of S epitopes associated with clonotypes was carried out by recognizing in the repertoires the TCRβ experimentally associated with S epitopes (public databases MIRA [20] and VDJdb [21]), as well as inferring the S-specificity of clonotypes using the GLIPH algorithm.
The pEptide tcR matchinG predictiOn (ERGO) tool [22] was used to classify TCR-peptide binding affinity toward S peptides.The CDR3 of TCR is the major determinant of T cell specificity.The complete list of viral peptides and CDR3 sequences used as input on ERGO is reported in Table 3, respectively.As output, ERGO produces a binding probability score (maximum value 1 if the TCR and the peptide bind and 0 otherwise).

Statistical analysis
Statistical tests were selected based on appropriate assumptions with respect to data distribution and variance characteristics; p values < 0.05 were considered statistically significant.Statistical significances are reported in the figure and/or the figure legend.Statistical tests were performed with GraphPad PRISM software 9.3 (GraphPad Software, La Jolla, CA, USA).

Characteristics of the enrolled patients and donors
We enrolled 14 consecutive patients tested positive for SARS-CoV-2 with mild to severe COVID-19 [23] from January 24th to July 7th, 2022.At that time the Omicron variants were predominant in our setting in Calabria region, Southern Italy [24].Patients with mild COVID-19 symptoms for a maximum of 7 days who did not require oxygen support presented to the center dedicated to early therapies for COVID-19 [25] to receive neutralizing monoclonal antibodies (moAbs) or antivirals (i.e., remdesivir, nirmatrelvir/ritonavir or molnupiravir).Patients with moderate or severe COVID-19 were admitted to hospital.Blood samples for the purpose of this study were collected at least 24 h before starting treatment with moAbs or antivirals.Among these patients, those who did not receive any doses of the approved vaccines were categorized as "not vaccinated" (NV), while those who received at least two doses of BNT162b2 mRNA COVID-19 vaccine were included in the "vaccinated" (V) group [26].Immunocompromised patients were defined as those affected by onco-hematological diseases, primary/ acquired immunodeficiency, systemic inflammatory diseases, or those who received rituximab, methotrexate and/or other immunosuppressive drugs.The characteristics of the participants are summarized in Table 1.

Diversity and clonality of COVID-19 TCRβ repertoires
We performed TCRβ sequencing of PBMCs samples both before (pre-stimulation TCRβ repertoires) and after (post-stimulation TCRβ repertoires) T cell enrichment following in vitro 12-day stimulation with a pool of S peptides (Fig. 1).Collectively, we obtained 66133 and 67160 distinct clonotypes (TCRβ with unique CDR3 amino acid sequence) from the pre-and post-stimulation TCRβ repertoires (Fig. 1).Information on TCRβ repertoires characteristics is reported in Table 2 and Supplemental data.
We first analyzed the diversity and clonality of prestimulation COVID-19 TCRβ repertoires, as compared to a control TCRβ dataset obtained from PBMCs sample of healthy individuals (n = 14) from the TCRB-V4b Control Database [27], that matched for age, gender, and ethnicity with our COVID-19 cohort.To rigorously evaluate the diversity of repertoires, we utilized TCRβ repertoire diversity and clonality estimators.Both the Gini-Simpson and D50 diversity indices were significantly lower in the COVID-19 groups compared to the control group, indicating less diverse COVID-19 repertoires (Fig. 2a).Conversely, the Gini coefficient was significantly higher in the COVID-19 groups than in the control group, indicating greater inequality of clonotype frequency (Fig. 2a).Overall, these TCRβ repertoire estimators demonstrated that COVID-19 TCRβ repertoires had higher clonality and lower diversity than the repertoires of healthy subjects, consistent with the expected clonal expansion typical of a cellular response to viral antigens.Notably, the two patients who developed severe COVID-19 exhibited very narrow true diversity and high clonality (Fig. 2b), consistent with previous studies [15].We also compared the distribution of Vβ gene segments between COVID-19 groups.In the V group, the TRBV7-2, TRB29-1, and TRBV30 Vβ gene segments showed a significantly higher clonal fraction than those observed in the NV (Fig. 2c).

In vitro expansion of S-specific T cells was equally efficient in PBMCs from vaccinated patients
The in vitro stimulation promotes the expansion of peptide-specific T cell populations, thus allowing the detection of low-frequent clones [20].To ensure that our S-specific T cell enrichment strategy was effective, we TCRβ repertoire analysis involved a filtering step to exclude confounding sequences not shared between pairs of pre-and post-repertoires.The mapping of S epitopes associated with clonotypes was carried out by recognizing in the repertoires the TCRβ experimentally associated with S epitopes (public databases MIRA [20] and VDJdb [21]), as well as inferring the S-specificity of clonotypes using the GLIPH algorithm [12] counted the T cell clones responsive to S-peptides using the ELISPOT assay (Fig. 3a, b).Compared to unstimulated PBMCs, those stimulated with S-peptides plus IL-2 showed a higher number of IFNγ-secreting T cell clones (Fig. 3a).The fraction of T cells responsive to S-peptide was not significantly different between V and NV individuals (Fig. 3b).
To further ascertain that post-stimulation TCR repertoires contained S-enriched clonotypes we took advantage of TCRβ sequences with known S-epitope specificity present in the public VDJdb [21] and MIRA [20] databases.The merged the VDJdb/MIRA dataset contained pairings of 26,422 unique TCRβ sequences and 388 S peptides from VDJdb and MIRA databases.The overlap of this dataset with our COVID-19 dataset resulted in 952 unique TCRβ sequences associated with 184 S-epitopes.For each patient, the cumulative frequency of S-specific TCRβ sequences was significantly increased in post-stimulation TCR repertoires (median 7.8; 95%CI 1.4-17.5)than pre-stimulation ones (median 0.8; 95%CI 0.2-1.9%,P < 0.0001, Wilcoxon matched-pairs signed rank test) (Table 2).Furthermore, the 10 most abundant clonotypes of each post-stimulation repertoire were increased in frequency compared to the pre-stimulation frequency (Fig. 3c, d).The S-specific TCRβ sequences accounted for an average of 12.9% (range, 0.7-38.8%)and 10.8% (range, 1.3-29.1%) of the entire post-stimulation repertoires of NV and V group, respectively (Table 2), indicating that the T-cell stimulation with S-peptides was similarly efficient in PBMCs from V and NV patients, in agreement with the ELISPOT results.

Epitope mapping of S-specific TCRβ clonotypes in COVID-19 repertoires
To further infer the specificity of the TCRβ clonotypes, we performed a clustering of TCRβ sequences based on sequence similarity to S-specific public clonotypes present in the VDJdb/MIRA dataset.To this end, the GLIPH2 algorithm can reliably group TCRs of common specificity from different T-cell samples, organizing clusters of TCR sequences according to their likely antigenic specificities [12,13,28].
We applied GLIPH2 to a dataset that included our COVID-19 TCRβ clonotypes and the VDJdb/MIRA dataset.We identified 347 specificity groups shared by N, NV, and VDJdb/MIRA datasets, while 21 and 33 specificity groups were shared by NV and VDJdb/MIRA, or V and VDJdb/MIRA datasets, respectively (Fig. 4a, b).The presence of TCRβ sequences from the VDJdb/MIRA dataset in each of the clustered specificity groups allowed us to infer the epitope specificity of the TCRβ sequences from the NV and V datasets within the same cluster [13] (Fig. 4c, d).We identified a total of 1171 TCRβ sequences associated with 149 peptides, of which 130 epitopes were shared between V or NV groups (Table 3).The cluster TCRβ sequences account for an average of 30% of the entire post-stimulation repertoires.Among the 130 shared peptides, TCRβ clonotypes associated with S regions S135-177, S264-276, S319-350 and S448-472 were significantly more abundant in the poststimulation repertoires of V group than the NV group (Fig. 5a, b and Table 4).Conversely, TCRβ clonotypes associated with protein region S673-699 was significantly more abundant in the post-stimulation repertoires of NV group than the V group (Fig. 5a, b and Table 4).Moreover, the S645-645 and S751-760 peptides were exclusively associated with TCRβ clonotypes in the NV group, while S778-789, S863-871 and S1260-1269 peptides were exclusively associated with TCRβ clonotypes in the V group (Fig. 5a, b and Table 4).These results indicate a distinct profile of TCR epitope specificity between the N and NV COVID-19 groups following infection with SARS-CoV-2.

Effect of epitope mutations on TCR-peptide binding affinity
The T cell expansion procedure was performed using ancestral Wuhan-Hu-1-derived S peptides to stimulate PBMCs from Omicron/BA2, BA.4, and BA.5-infected patients.One potential concern was that TCR specificity might vary between individuals exposed to both the Wuhan-Hu-1 and Omicron antigens (group V) versus those exposed only to the Omicron antigen (group NV).
Among the relevant S peptides of the ancestral Wuhan-Hu-1 strain, the S673-688 and S975-985 were subjected to amino acid substitution within the variants BA.1, BA.2, BA.2.12.1, BA.2.75, BA.4,BA.5, BQ1.1, XBB1.5.We addressed the impact of amino acid substitution on TCRβ binding by performing an in-silico analysis using ERGO, a reliable and robust tool for TCR-peptide binding affinity prediction and classification [22].ERGO-II results indicated that the variant in the S975-985 region (S981L > F, of Omicron/BA.1) increased the peptide binding probability score of TCRβ clonotypes from both N and NV groups, while the variant in the S673-688 region (S679N > K, S682N > K of Omicron/BA2, BA.4,BA.5, BQ1.1 variants) did not affect the peptide binding probability score of TCRβ clonotypes associated with this region (Fig. 5c).These results suggested that epitope mutations did not preclude cross-reactive recognition by TCRβ clonotypes of our groups.

Discussion
In this study, we first highlighted the differences in the general characteristics of the TCRβ repertoire between individuals with COVID-19 who were either vaccinated or unvaccinated.Then, we specifically examined the S-specific TCRβ repertoire, identifying differences that may be important in predicting disease progression and cross-reactivity towards viral variants.
The ability of adaptive immunity to achieve effective TCR diversification may determine the chances of improved outcomes and immune control in infectious and cancer diseases [29][30][31][32].A study on the T cell repertoire in SARS-CoV-2 infected patients showed that TCR diversity may influence disease outcomes, and that patients with asymptomatic or mild clinical infection have a highly diversified TCR repertoire, while patients with severe COVID-19 have a less diverse TCR repertoire [15].Consistent with this previous observation, we found that patients with SARS-CoV-2 exhibited lower diversity and higher clonality compared to the healthy control group, which is expected in the presence of a productive infection.Furthermore, COVID-19 patients with breakthrough infections after vaccination, all of whom had a mild clinical course, showed a trend towards greater diversity compared to non-vaccinated COVID-19 patients, although this difference did not reach statistical significance.Additionally, the distribution of diversity indices within the two groups appeared different, with much more homogeneous diversity values in vaccinated patients compared to the unvaccinated ones.It is possible that in patients with breakthrough infections, the presence of established immunological memory from vaccination allowed for a more expansion of S-specific T cells, resulting in a repertoire with a greater and homogeneous diversity compared to the unvaccinated group.It should be noted that the two patients who developed severe COVID-19 exhibited very narrow true diversity and high clonality, consistent with the previous referenced studies [15].
The analysis of Vβ gene segments distribution within the TCRβ repertoire showed that TRBV7-2, TRB29-1, and TRBV30 Vβ gene segments had a significantly higher clonal fraction than those observed in the unvaccinated group.T cells expressing the TRBV7-2 gene segment have been associated with various clinical conditions, including immune response to human rhinovirus infection and autoimmune diseases such as multiple sclerosis and rheumatoid arthritis [33,34].Our results further expand the conditions associated with TRBV7-2 Vβ gene segment and support further investigations to uncover the clinical significance of this association in COVID-19.
The S673-699 region was specifically linked to TCRβ clusters in the unvaccinated group, among the S regions that were differentially associated with TCRβ from vaccinated or unvaccinated groups.Some researchers have suggested that this region may have super antigenic properties, which could potentially contribute to the severe immune response observed in some COVID-19 patients, leading to cytokine storms and multiorgan failure [35][36][37].Superantigen-mediated T cell expansion occurs through a less specific interaction with TCRs, leading to the activation of a large proportion of the T-cell population [38].Although the super antigenic character of the S673-699 region is still being investigated, it is possible that the efficacy of the mRNA COVID-19 vaccines in preventing severe disease and death is also due to reduced exposure to the super antigenic viral    determinant, through neutralizing antibodies or reduced viral replication.We also found that TCRβ clonotypes associated with S regions S135-177, S264-276, S319-350 and S448-472 were significantly more abundant in the vaccinated group than in the unvaccinated group.There could be several reasons for this difference, which are likely related to the various ways in which these peptide regions are expressed and presented to the immune system between vaccination and natural infection.In our small cohort study, the presence of specific clonotypes associated with these S regions was found to be linked with a milder course of the disease.Therefore, detecting T cells that are specific to these regions may have a positive impact on disease prognosis and could potentially be used as a prognostic tool.
Our study obtained TCRβ sequences specific for the S protein from PBMCs of patients infected with the Omicron/BA2, BA.4, and BA.5 variants, after stimulation  3) were evaluated in-silico for predicted binding affinity by the ERGO tool (Mann Whitney test) with a pool of S peptides derived from the ancestral Wuhan-Hu-1 sequence.In vitro expansion of these cells indicated that TCRβ cross-reacted with peptides from the ancestral strain.Our in-silico analysis using the ERGO tool showed that the variant in S975-985 region increased the probability score of TCRβ clonotype peptide binding in both vaccinated and unvaccinated groups, while the variant in the S673-696 region did not affect the probability score of peptide binding of TCRβ clonotypes associated with this region.These findings suggest that TCRβ clonotypes from both groups can recognize crossreactive epitopes despite mutations.
The study provides TCRβ sequences, which is a valuable information for understanding the immune response.However, it is important to note that the lack of alpha chain information represents a limitation to the study.Indeed, the alpha chain plays a critical role in shaping the T cell receptor's specificity and affinity for antigen recognition.However, although the absence of alpha chain information may hinder the interpretation of epitope specificity, this limitation is diminished by the fact that the prediction algorithms we used were predominantly trained with TCRβ data.Moreover, the potential confounding factors in our study, particularly in relation to the small and diverse sample size, limit the ability to draw any association of TCR signature with clinical variables.The presence of both immunocompetent and immunodeficient patients, the limited number of severe cases (n = 2), and variations in vaccination status and HLA types indeed pose challenges in interpreting the TCR data with high precision.Moreover, we recognize that the limited scope of our sample size, comprising only 14 individuals, restricts our ability to generalize our findings to a broader population.This sample size limitation is especially pertinent given the complexity of TCR repertoires and their interactions with various clinical variables.Our study's focus on selected aspects of the immune response further narrows the scope of our findings.

Fig. 1
Fig.1Study workflow.The PBMCs from a whole blood venous sample of COVID-19 patients were divided into two aliquots.One aliquot was directly subjected to TCRβ sequencing (pre-stimulation repertoires).The other aliquot was first stimulated with S-peptides pool, in the presence of IL-2 for 12 days, to promote the expansion of S-specific cells, and then was also subjected to TCRβ sequencing (post-stimulation repertoires).TCRβ repertoire analysis involved a filtering step to exclude confounding sequences not shared between pairs of pre-and post-repertoires.The mapping of S epitopes associated with clonotypes was carried out by recognizing in the repertoires the TCRβ experimentally associated with S epitopes (public databases MIRA[20] and VDJdb[21]), as well as inferring the S-specificity of clonotypes using the GLIPH algorithm[12]

Fig. 2
Fig. 2 TCRβ repertoires diversity and clonality of COVID-19 and healthy groups.a For Gini-Simpson, D50 and Gini diversity significance see the Methods section.Statistical comparison was performed by Kruskal-Wallis test.b TCR repertoire fraction of the top 10, rare, and hyper-expanded clonotypes for all repertoires.Top 10 clonotype refers to the 10 most abundant clonotypes from each repertoire; rare and hyper-expanded refer to clonotype with a frequency less than 10 -5 , or greater than 10 -2 , respectively.The position of the value observed for the patient who experienced severe COVID-19 is highlighted in yellow.c Comparison of Vβ usage between V and NV COVID-19 groups.Statistical comparison was performed by Mann Whitney test

Fig. 3
Fig. 3 In vitro expansion of S-specific T cells.a Representative IFNγ-ELISPOT results for a case of NV and V COVID-19 patient.b Frequency of SFU of IFNγ-secreting cells following stimulation with the S-peptide pool in NV and V COVID-19 patient.*p = 0.003, **p = 0.006, Fisher's test).Statistical comparison was performed by Kruskal-Wallis test.c, d Representative examples of tracking of the top 10 most abundant clonotypes from the post-stimulation repertoires of NV (c) and V (d) COVID-19 patients

Fig. 4
Fig. 4 Inference of TCRβ clonotypes specificity by GLIPH analysis (a) Network analysis of the most significant (Fisher score < 10 -4 ) specificity groups annotated with TCRβ sequences from NV (red), V (green), and VDJdb/MIRA (blue) datasets.Each dot is a specificity group, edges indicate the presence of identical TCRβ sequence(s) shared across two specificity groups.b Shared number of S-epitopes in VDJdb/MIRA, V and NV groups.c, d Characteristics of a representative clustered specificity group, including a representation of amino acid sequence alignment of both the multiple TCRβ in the cluster (c), and the associated S region (d)

Fig. 5
Fig. 5 Mapping of TCRβ associated with S epitopes as resulted from GLIPH2 analysis.a Median clonal fraction refers to the median of clonotype frequencies from the individual pre-stimulation repertoires.Dimension of bubbles correlates with the number of distinct clonotypes associated with the epitope.b Glimpse of panel A highlighting the S regions (overlapping epitopes) associated to TCRβ differentially abundant between COVID-19 groups.C TCR binding prediction of S peptides from original Wuhan-Hu-1 strain and VOCs variants.The distinct TCRβ sequences associated with the indicated S peptide (Table3) were evaluated in-silico for predicted binding affinity by the ERGO tool (Mann Whitney test)

Table 4 8 c
Summary of TCRβ associated with spike peptides.In bold, the data represented in Fig. 5b are highlighted a Amino acid position in the spike protein (uniprot ID P0DTC2) b Number of clonotypes associated to epitope within the indicated spike region; showing a post/pre frequency > Median of the enrichment (post/pre stimulation ratio of clonal fration) showed by the clonotypes associated to epitope within the indicated spike region d Median of the pre-stimulation clonal fration showed by the clonotypes associated to epitope within the indicated spike region e Statistical differences were evaluated by Mann-Whitney test

Table 2
Information on immune repertoire sequencing and analysis a n_cln-number of unique CDR3β amino acid sequences b Number of clonotypes shared with the MIRA/Adaptive dataset c Fraction of repertoire' sequences relative to S-specific clonotypes

Table 3
Characteristics of the specificity clusters obtained by GLIPH analysis

Table 3
(continued) Fisher_score estimates the bias of the pattern presenting in target data versus reference data using Fisher exact test; b Patter of the cluster; global pattern contains '%' , which indicates position allowing variants; local pattern starts with 'motifc Number of distinct CDR3s included in the cluster; d Fisher_score estimates the bias of the pattern presenting in target data versus reference data using Fisher exact test [12]; e CDR3 amino acid sequences of clonotypes from the COVID-19 dataset; the clonotypes from VDJdb/MIRA dataset; f S-peptides associated with clonotypes of the VDJdb/MIRA dataset included in the cluster; g Amino acid position in the spike protein (uniprot ID P0DTC2)